LIEB-THIRRING INEQUALITY FOR A MODEL OF PARTICLES 
WITH POINT INTERACTIONS 



RUPERT L. FRANK AND ROBERT SEIRINGER 

Abstract. We consider a model of quantum-mechanical particles interacting via 
point interactions of infinite scattering length. In the case of fermions we prove a 
Lieb-Thirring inequality for the energy, i.e., we show that the energy is bounded 
from below by a constant times the integral of the particle density to the power |. 
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1. Introduction 

Quantum-mechanical models of particles with point interactions have become rel- 
evant in cold-atom physics, where one encounters systems where the range of the 
interactions among the atoms can be much shorter than their average distance, while 
the scattering length can be much larger. In the limit of zero interaction range and 
infinite scattering length, which is referred to as the "unitary limit" , one is left with a 
system without intrinsic length scale. For further discussion of this topic we refer to 
[HI El m ESI n] and, in particular, to the articles in [2i] . 

While the two-body problem for particles interacting via point interactions is well 
understood [l], it remains an open problem to establish the existence of a model for 
N > 2 particles with only two-body point interactions (see [H [10] and references 
there). Such a model can only be expected to exist for spin ^ fermions, due to the 
Efimov effect [9j, i.e., the existence of three-body bound states despite the absence 
of two-body bound states. We consider here a model with more complicated point 
interactions which is manifestly well-defined, in the sense the quadratic form for the 
energy is positive. This model is well-defined even for bosons, and does not allow for 
any bound states even in this case. 

We shall prove that the model we consider satisfies a Lieb-Thirring inequality, i.e., 
the energy of fermions is bounded from below a semiclassical approximation to the 
kinetic energy. Up to the value of the constant, this inequality is the same as the one 
obtained by Lieb and Thirring |19^ I20j for the kinetic energy only, i.e., in the absence 
of any interaction. We recall that this inequality of Lieb and Thirring played a crucial 
role in a short and elegant proof of the stability of matter [17], which had been proved 
earlier by Dyson and Lenard [8] by different means. 
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Since our model concerns interacting particles, we cannot use the method of Lieb 
and Thirring to reduce the problem to the spectral analysis of a one-body operator. 
Instead, our strategy is closer in spirit to the work of Dyson and Lenard [8j, where 
the Pauli principle enters only via a local exclusion principle, which implies that the 
local kinetic energy cannot be zero for more than one particle. A similar strategy was 
recently employed in [2T] to study a system of anyons in two dimensions, and ideas 
from [21] are crucial to our proof. 



2. Model and Main Result 
Let iV > 2, X = (xi, . . . , xat) G M^^ and let g : M?^ M denote the function 

l<i<j<N ' ' ■'' 
Our model is defined via the quadratic form 

N 

Qif) = Y. 9(.Xf\Vif{X)\'dX (2) 

on L'^{M.^^ , g{X)'^dX), where Vj stands for the gradient with respect to Xi G M^, and 
dX = Uk=idxk- The norm on the space L'^{M.^^ , g{X)'^dX) will simply be denoted 
by II • II . The model ([2]) was introduced in [2|, where it was shown that the quadratic 
form Q{f) gives rise to a non- negative self-adjoint operator with purely absolutely 
continuous spectrum [0,cxd). 

We denote by C L'^{M.^^ , g{X)'^dX) those functions / that have the property 
that there is a partition of {!,... , N} into q disjoint subsets such that / is antisym- 
metric in the variables corresponding to each subset. In particular, Ai are the totally 
antisymmetric functions, while A^ are ah functions in L'^{R^^ , g{X)'^dX). Note that 
A^ C -^^1- The functions / G A^ thus represent the spatial part of totally anti- 
symmetric functions of space and an internal degree of freedom ("spin"), where the 
spin is allowed to take q different values. 

For given / G L'^{R^^ , g{XydX) we define its density gj G L^{R^) by 

TV 

^?/(^) = E/ 9{X,f\f{X.)\'dX, (3) 

with Xi = (xi, . . . , Xj_i, X, Xj+i, . . . , xjy) and dXi = Hj^i dxj. Our main result is the 
following Lieb-Thirring type inequality. 

THEOREM 1. For some constant C > (independent of N and q) we have 

Q{f) >4rJ Bfixf'dx (4) 

for all f e H\R^^,g{XfdX)nA^ with \\f\\ = 1. 

As mentioned in the introduction, Lieb and Thirring proved ^ in the case g ^ 1. 
Improved bounds on the constant C in this case were obtained in [7] (see also [El [13]). 
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For C open and with finite measure, the ground state energy for N particles 
confined to O equals 

E{N,q,n)=M{Q{f) : f e C^in'') H , ||/|| = l} (5) 
for our model. An immediate corollary of ^ and Holder's inequality is that 

C A^5/3 

where |r2| denotes the volume of il. 

Remark 1. With A^, denoting the Laplacian with respect to the variables x^, we have 
Afe g{X) = for all X E M^iv ^^^j^ |^. _ ^.^ > g for all 1 < f < j < A^. If we 

define 

T, = {X€ R^^ : \xi -Xj\>eyi<i<j < N} (7) 

for e > 0, then an integration by parts gives 

Af » ^ r 

g{xf\V^f{xfdX = Y, \V^g{X)f{xfdX 

i=l "^"^E i=l "^"^s 



(8) 



i<i< 

where dS* denotes the induced surface measure on dT^ , the boundary of . Hence our 
model indeed describes particles with point interactions supported on the union of the 
hyperplanes defined by Xi = xj for 1 < i < j < N . 

Note that the last term in ([8]) is negative and vanishes in the limit e — )• if |/(A)p 
vanishes faster then linearly on the hyperplanes where Xi = Xj. In particular, if 
f{X) = ^(X)/(/(A) for ^' G C^{R^^), then 

N 

In general, the functions / need not vanish for coinciding arguments, however, and 
the energy can be lowered due to the attractive nature of the last term in ([8]). 

Remark 2. If we replace l/|x| in ([T]) by a smooth, strictly positive functions (p{x) 
and define, accordingly, = Ylii<i<j<N ^i^i ^^"^ ^(^) = 9{^)f{^)i then an 
integration by parts gives 

N N 

Y g{X?\V,f{X)\^dX = Y \V^^{XfdX 

+ [ |vI,(X)|2££iMWdA. (10) 
9{X) 



The effective potential g{X) ^ 'Ylik=i ^kg{X.) can alternatively be written as 

^ 9{X) ^ ' 

l<i<j<N ' 
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and thus contains not only two-body terms but terms involving arbitrarily many par- 
ticles. In particular, the presence of other particles besides i and j decreases the 
summands in (jlip and hence weakens the interaction. This weakening leads to im- 
proved stability properties as compared to the case with only two-body interaction. 



We note that our main result Theorem [T] holds for more general functions g{X). As 
an example, one could replace g{X) by 

■ _l 1' 



E 

l<i<j<iV 



(12) 



for a < 0, corresponding to point interactions with finite scattering length a. A Lieb- 
Thirring inequality holds also in this case, with a constant C that is bounded from 
below uniformly in a. As a — )• one obtains non-interacting particles. 

Our method can also be applied to analogous models in dimensions n > 3, where 
the function l/\x\ in ([T]) should be replaced by the Green's function of the 

Laplacian in M*^. 

The remainder of this paper is concerned with the proof of Theorem [TJ The proof 
will be split into two main parts. In the first part in Section [3l we prove a local 
exclusion principle (Prop. [1]) which states that the energy of n particles in a finite 
cube is strictly positive for n > q and grows at least like n — q, in fact. The second 
part in Section [J] concerns the proof that the energy Q{f) dominates the L^(R^) norm 
of Vy/gJ. The inequality we prove in Prop. [2] is, up to a constant, equal to a well- 
known inequality by the Hoffmann-Ostenhofs [12] in the case without interaction. In 
Section [5] we demonstrate how to obtain Theorem [1] from Propositions [1] and [21 The 
strategy of the proof is similar to the one in |21] where a Lieb-Thirring inequality was 
proved for anyons in two dimensions. 

3. Local Exclusion Principle 

In this section we shall prove the following. 

Proposition 1. Let f £ H\R-^^ , g{X)^dX) n with \\f\\ = 1, and let d C be 
a cube of side length L > 0. Then 

^^^^<7(^)2|V,/(X)|2;^cJx,)dX> Aj^jT gf{x)dx-q^ (13) 

for a constant k > independent of N , f , L and the location of the cube Cl. 

Here and in the following, XQ denotes the characteristic function of a set Q C M'^. 
The constant k appearing in (I13p is the same as the one in Lemma [2] below. 
We note that this proposition implies that, for any integer M > 1, 

NM"^ - qM^ 

E{N,q,CL)>k'^^^!-^j^. (14) 

To prove ([H]), simply divide the cube Cl into disjoint cubes with side length 
L/M and apply Proposition [1] for every cube. The optimal choice of M is close to 
{5N/2q)^/^ . In particular, the maximum over M of the right side of (|14p is proportional 
to kN^''^/{L'^q^l^) for large N. 
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The proof of Proposition [T] will be divided into several lemmas. The following lemma 
is well-known [5j; we include its proof for the convenience of the reader. 

Lemma 1. Let Q be an open and convex subset o/M" with diameter d < oo, and let u 
be a strictly positive function in L^{Q). For all f G H^{Q,,uj{x)dx) with f{x)dx = Q, 
we have 

|V/(x)|2a;(x)d2; > ^— ,^ ' ' , / \f{x)\'^io(x)dx, (15) 

where \W^\ denotes the volume of the unit ball in M" and 

= \nr^n ( \\ ~f r\ (16) 

with Bf.{x) denoting the ball of radius r centered at x. 
Proof. Since O is convex, we can write 

f{x) - f{y) = [\x- y)Vf{tx + (1 - t)y)dt . (17) 







Integrating this identity over y G and changing variables to z = + (1 — t)y gives 
Since \x — z\ = {1 — t)\x — y| < (1 — t)d for x,y £ ^l, we obtain the bound 

1 r rl-\x-z\/d 1 

< ^ / \x-z\Nf(z)\ / -—dtdz 

IM )\- y^l MM )\ (l-t)n+l 



< 



f 1 



n\^\ Jq \x — z\ 

For X, z € we have 



-\Vf{z)\dz. (19) 



|^_l|n-i ^XBAx)iz)dr + ^ . (20) 

Hence ()19p can be rewritten as 



{n-l) f \ f \Vf{z)\dzdr+-^ [ \Vfiz)\dz 

Jo ^ JnnBr(x) " Jq 



We introduce the maximal function 



(21) 



, , W(x) |V/(^)|c^(^)rf^ 
ma;(x) = sup -— . (22) 

Then 

\Vf{z)\dz < \M''\r''M^m^{x) (23) 



innBr(x) 

with Mu defined in In particular, from ([21]) we obtain the bound 

jn+l 

|/(x)| < ^|B"|M,m,(x), (24) 
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and thus 

J2(n+1) j- 

\f{x)\'^u{x)dx<——\WfMl / m^{xfuj{x)dx. (25) 
We further define 

Nuj = sup — ' ' , (26) 

r>o,xen jQnBrix) ^\z)az 

and note that < 3"Ai^. Proceeding as in the proof of Theorem 1(c) in Sect. 1.3.1 
of [22| we see that 



m^{xyuj{x)dx <2^N^ \V f{x)\^uj{x)dx . (27) 
! Jn 

In combination with (I25p this proves our claim. □ 

Let now C = (0, 1)^ C denote the unit cube in M^. 

Lemma 2. Assume that f G i7-'^(M'^xM'^) is antisymmetric, i.e., f{xi,X2) = —f{x2,xi) 
for xi,X2 G K'^. For almost every (x^, . . . ,xn) G we have 

2 

V/ g{Xf\Vif{xi,X2)\^dxidx2>2k[ g{Xf\f{xuX2)\^dxidx2 (28) 

for a constant k > that does not depend on N or the variables {x^, . . . ,xn)- 

Proof. Since / is antisymmetric, it has zero average with respect to the variables 
(xi, X2). Hence we can apply the previous lemma, with = C xC C and a;(xi, X2) = 
g(X)^. The claim is thus proved if we can show that the corresponding in ()16p is 
bounded independently of and (X3, . . . , x^r). 

Consider a ball S C of radius r centered at some point (wi, ^2) G C xC. Without 
loss of generality, we may assume that r < ^/6, for otherwise C x C C B and hence 
the expression after the sup in (I16p is independent of r, and is thus the same as for 
r = V6. We can get a lower bound on g{X) by taking the minimum in each term in 
the sum in ([1]). This gives 

Iqiu — 7;jo -1- 9r ^ — ^ \ — T'j -i- T 7/I0 — T'j -4- r I ^ — ^ 



wi — W2 \ + 2r ^ \ \ w\ — xA+ r \w2 — x, + r / ^ x,- — x; 



(29) 



for (xi,X2) G i?. On the other hand, simple estimates yield 



1 r'° 

-dx\dx2 < const j ^ (30) 



b|xi-X2P (Itfi - t(;2| + 2r)2 

and 

f I r^ 

/ 1 n 7dxidx2 < const j — ^ j r- (31) 

Jb \xi - Xj\\xi - xi\ {\wi -Xj\+ 2r)[\wi - x/| + 2r) 

for j, I > 3, etc. These imply, in particular, that 

/ 5(X)2dxidx2 < const inf g{X) . (32) 

Jb (xi,x2)eB 

Since is a cube in M^, we also have |r2 n S| > 2~®r® for r < 1, and hence is 
bounded independently of and (xs, . . . , xat). □ 
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Lemma 3. Let A C {!,..., A^}, Xa = {xijigA o-nd dXA = H/eA^^'- -^'^'^ / ^ 
H'^{^^^ ,g{XfdX) n we have, for almost every X \ Xa, 



V / g{Xf \V,f{X f dXA > k {\A\ -q) [ g{Xf\f{X)\''dXA 



(33) 



with k as in Lemma\^ 



Proof. The set A can be divided into q subsets Bj (some of which may be empty) such 
that / is antisymmetric in the variables in each subset. We may assume that |yl| > q, 
in which case at least one such subset Bj contains more than one element. We shall 
show that for all such Bj 

!^Q{Xf\V,f{X)\'dXB,>k{\B,\-l) [ g{Xf\f{X)\^dXB^ (34) 

for each fixed X \ Xb^ Summing over j this implies the result. 

Inequality (|34p follows immediately from Lemma [51 noting that the sum over \Bj\ = 
n>2 elements can be written as l/[2(n— 1)] times a sum over all ordered pairs, yielding 
the inequality with \Bj \ — 1 replaced by |-Bj|. □ 

Since g{X) is invariant under translation and rotation of all the coordinates Xj, 
inequality (j33p clearly holds for all unit cubes in R'^, irrespective of their location or 
orientation. A simple scaling argument shows that it actually holds for all cubes of 
side length L > if the right side is divided by L?' . 

Proof of Proposition [1]. Using that 

1= E \{xc,{xi)J{{l-XcA^i)) (35) 
A(Z{l,...,N} leA l^A 

we have 

N 

/ g{Xf\V,f{X)\'xc,{xi)dX 

= E E /.,^?wi^^/wi'n^c.(^')n(i-xc,(xo)dx. (36) 

Ac{l,...,iV} ieA •^'^^^ l£A 1<^A 

Applying the previous lemma (and the remark after its proof), we obtain the lower 
bound 

A Y m-q)j g{Xf\f{X)\^\{xcM)\{{^-XcM))dX, (37) 
Ac{l,...,Af} leA l^A 

which is equal to the right side of (|13p . □ 



4. Density Bound 

Recall that the one-particle density Qf of a function / G L'^(M.'^^ , g{X)'^dX) is 
defined in ([3]). 
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Proposition 2. Let f £ H^{R^^ , g{XfdX). Then its one-particle density Qf satisfies 
^ £ H'^{R^,dx) and 

2 



Qif) > Q / Vj£./(x) dx. (38) 
It is well known that in case g = 1, ()38|) holds with the prefactor ^ replaced by 1 

ni. 



Proof. We first show that it is enough to prove the inequality for symmetric functions 
/. For any /, define / by 

1/2 

I ^ 1. I , (39) 



where Sn denotes the permutation group. The function / is obviously symmetric, and 
Qf = Qj. It is easy to see that the map |/p i— )• Q{\f\) is convex (see, e.g., Lemma A.l 
in [18]), hence 



Q(/)<g(|/|)<Q(/). (40) 

In particular, it is enough to prove ([38|) for symmetric /. 

Assume now that / is symmetric. Then the density Qf can alternatively be written 

as 



Qf{xi) = N / g{XY\f{X)\'dx2 ■■■dxN 

Jr3{N-1) 

= N{N - 1) / ^-—\f{X)\^dX2 ■■■dXN 

Jr3(]V-i) \Xi - 

+ N{N -l){N-2) [ -\f{X)\^dx2 ■■■dxN 

JkHN-1) \Xi - X2\ \xi - xsl 

+ 2N{N -1){N -2) [ -\f{X)\'^dx2---dxN 

7k3{JV-i) \Xi - X2\ \X2 - Xsl 

+ iV(iV- l)(iV-2)(iV-3) / r-> 7\f{X)\'^dx2---dxN 

7R3(iV-i) \XI - X2\ \X3 - X4| 

+ N [ R{x2,...,XN)\f{X)\^dx2---dxN (41) 

Jr3{JV-1) 

where R denotes the sum of all the terms in g not depending on xi. Let us denote 
the various integrals on the right side by Qi, ■ ■ ■ , Q5, and the prefactors in front of the 
integrals by ni, . . . , 725. That is, Qf = X]j=i "-j^'j- 

Since l/|xi — X2I is invariant under translations of both xi and X2, we have 

^|/(X)|2 = 2Re^ 

\Xl — X2\ \Xi — X2\ 

Integrating this identity over X2, . . . , Xjv, we obtain 

ViQi{xi) = 2Re [ -. ^-—f{xy{Vi + V2)f{X)dx2---dxN. (43) 

J |xi-x2r 

The Schwarz inequality then implies that 

|Vi£>i(xi)|2 <4£>i(xi) / L^|(Vi + V2)/(X)|2dx2---dx^. (44) 

J X1-X2K 



(Vi + "^2)72 ir^l/ Wl' = 2Re I _ .^ f{XnVi + V2)f{X) . (42) 
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Jr3 4 J^3 Qi{x) 

<[ I ^|(Vi + V2)/pdX<4 / L^|Vi/|2dX, (45) 

J^3N \Xl - X2Y Jr3N \Xl - X2\^ 

where we have again used the symmetry of / in the last step. 

In the same way one proceeds for the remaining parts of the density. The result is 
that 

/ |VV^|2< [ ^|(Vi + V2 + V3)/|'(iX 

7r3 Jr3n \xi - X2I \xi - xsl 

< 3 /" 1 ri r|Vi/|^dX 

J^3N \Xl - X2I \Xl - X-i\ 

+ 6/ ^ r|V2/pdX, (46) 

J^3N \Xl - X2\ \Xl - X3\ 
[ |V7^|2< / ^|(Vi + V2 + V3)/|'dX 

Jr3 JlR3iV \Xl - X2I \X2 - X3I 

< 6 /" 1 — ri — rlVi/pdX 

J^3N Xl - X2 X2 - X3 



+ 3/ -. i ^ r|V2/|2dX, (47) 

JrSN \Xi - X2I |X2 - X3I 

[ |VV^|2< / -. r-. A{Vl + V2)f\^dX 

Jr3 jR3iV |Xi - X2I |X3 - X4I 

<^f 1 — ri — - — r|Vi/pd^ (48) 

Jr3JV |Xi - X2I |X3 - X4I 



and 



[ |VV^|2</ R{x2,...,XN)\Vif\^dX. 



(49) 



Summing up, using again the symmetry of / and convexity of the map g J |Vy^p, 
this gives 

5 

/ \y^Q?<y,nJ \VM'<m[ g{Xf\VJ\^dX = 9Q{f). (50) 

Jr« ~[ iR3 ■' jR3iV 

This completes the proof of our claim. □ 

5. Proof of Theorem [T] 

We will deduce Theorem [1] from Propositions [1] and [21 following a similar strategy 
as in |21] . 

First, we note that it is enough to prove the theorem for N > 2q. li N < 2q, 
the statement can be easily deduced from Proposition [2] alone. In fact, applying the 
Sobolev inequality IIV99II2 > 'S'llv^Hg and Holder's inequality, we have 

Qif) > liNvojwi > f ii^/iis > f > f 9f{-r'dx (51) 
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in this case. 

Assume now that N > 2q. For given / G 
can choose a cube Qo C M'^ such that 



g{XydX) n A'^ and e > 0, we 

/ Qf{x)dx>2q and / Qf{xf'^dx>{l-e) Qf{xf'^dx. (52) 
Jqo Jqo J^^ 

Suppose that Qq is divided into finitely many disjoint cubes Qj. With the aid of 

Proposition [H we can bound 



«/) > EE/.„ ^(^f iv,/(x)i^x«.(..M^ ^-^lowii 



i j=l 



Qf{x)dx 



(53) 

where = max{t,0} denotes the positive part of a number t G M. Note that each 
summand on the left side is obviously non-negative, hence we can use the positive part 
on the right side. Moreover, from Proposition [21 we obtain the bound 

'^dx. (54) 



/ 


V^Qfix) 







In combination, (|53p and ()54p imply that 
Xk 



Qif) > E 



Qf{x)da 



+ 



1-A 



9 



VJgfix) 



dx 



(55) 



for any < A < 1. 

To construct the division of Qq, we proceed as follows [21]. Divide the cube Qo into 
8 disjoint cubes of half the size. If the integral of gj over one of the subcubes is less 
than 2q, the subcube will not be divided further and will be marked A. If all subcubes 
are marked A, then the division is undone and the cube Qo is marked B. For all the 
subcubes with integral Qf bigger or equal to 2q, we iterate this procedure. At the 
end, Qo is thus covered by finitely many disjoint subcubes of type either A or B. 

On every subcube Qi marked B, we have 2q < Jq gf{x)dx < 16q. This implies, in 
particular, that 

1 



gf{x)dx — q 



> 



gf{x)dx . 



The Poincare-Sobolev inequality on a cube [161 Thm. 8.12] yields 



VJgf{x) 



dx > S 



Qf - \Qi 



Qi 



(56) 



(57) 



(58) 



for some constant S* > independent of the size or location of Qi. The triangle 
inequality in L^{Qi) and a simple Schwarz inequality further imply that 

>^ii^/ii3-iQr^/^/ 

L6{Q,) ^ JQ^ 

An application of Holder's inequality finally gives 

5/3 



^fQj-\Qi\ ^ / 

JQr 



VJgjix] 



dx > 



r 5/ 



Q 



\ 2/3 



9f 



(59) 
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If we choose < A < 1 in such a way that 

i.e., A = (1 + 9k/{2S))~^ , we conclude that the contribution of a i? cube to the energy 
in (f55l) is bigger or equal to 

2-11/3 I 



k ^ S ^ 



Qf{xf/^dx. (61) 



Consider now a cube labeled A, where Jq. Qf < 2g. Pick a k > 2 and assume, for 
the moment, that 

/ 9f>^m-'/'([ QfY' ■ (62) 
In this case, it follows from (j59p that 

r 5/3 

/^Jvy^f..>|,«-2)^. (03) 

Hence, for our choice of A, the contribution from such an A cube to the energy in ()55p 
is at least 



We are left with studying those A cubes where 

5/3 



Jq, \Jq, 



(65) 



5/3 

We shall show that their contribution to the total integral of Qj^ is dominated by the 
contribution of all the B cubes. In order to see this, we note that our decomposition 
of Qo into subcubes can be organized in a tree. (Compare with Fig. 3 in [H].) Every 
A cube can be associated with a B cube, namely if it can be found by going back in 
the tree from the B cube, possibly all the way to Qo; Sind then one step forward. This 
allows to divide the A cubes into groups labeled by the B cubes. This division is not 
unique, in general, but this is not important; the only thing that matters is that every 
A cube can be associated with a B cube in this way. Pick a B cube (call it Qb) at 
level / G N of the tree, and let A{Qb) be the set of all those associated A cubes that 
satisfy (I65p . Since at every level 1 < j < / of the tree there are at most 7 A cubes in 
•^{Qb), we can bound 

On the other hand, for the associated B cube Qb have 



5/3 



\}Q,ei) ^ (2g)5/3 _ (2g)5/3 
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and thus 



Qi&A{QB) 



B 



In particular, the contribution to J q^^^^ of all the A cubes satisfying ([65]) is bounded 
by 28k/3 times the contribution of ah the B cubes. In other words, the contribution to 
/ Q^j^ of all the B cubes is bounded from below by (1+28k/3)-^ times the contribution 
of both the B cubes and the A cubes satisfying ([651) . 

After summing over all cubes, we thus get the lower bound 

Q{f)>-^[ 0f{xf'dx>-^{l-e)[ Qfixf'dx, (69) 

q^'-^ Joa q^'-^ Jr3 



with 



-11/3 . f 1 / 2\1 2-2/3 239-^/56977 



^ = ^^supmin<^— l--n>=^^ >0. (70) 



Since e > was arbitrary, this proves (jH with a constant 



C = min|c,^^|. (71) 
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